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Abstract. In recent work, the second author and various collaborators have shown 
using Evans function/refined semigroup techniques that, under very general circum- 
stances, the problems of determining one- or multi-dimensional nonlinear stability of 
a smooth shock profile may be reduced to that of determining spectral stability of the 
corresponding linearized operator about the wave. It is expected that this condition 
should in general be analytically verifiable in the case of small amplitude profiles, but 
this has so far been shown only on a case-by-case basis using clever (and difficult to 
generalize) energy estimates. Here, we describe how the same set of Evans function 
tools that were used to accomplish the original reduction can be used to show also 
small-amplitude spectral stability by a direct and readily generalizable procedure. 
This approach both recovers the results obtained by energy methods, and yields new 
results not previously obtainable. In particular, we establish one-dimensional sta- 
bility of small amplitude relaxation profiles, completing the Evans function program 
set out in Mascia&Zumbrun [MZ.l]. Multidimensional stability of small amplitude 
viscous profiles will be addressed in a companion paper [PZ], completing the program 
of Zumbrun [Z.3]. 



Section 1. Introduction 

In this paper, we study one-dimensional spectral stability in the small amplitude 
limit of smooth shock profiles, i.e., traveling wave solutions 

(1.1) u = u(x — st), lim u(z) = u±, 

z — >±oo 

arising under various regularizations of a hyperbolic system of conservation laws 

(1.2) u t + f(u) x = 0: 
specifically, viscous conservation laws of form 

(1-3) u t + f{u) x = (B(u)u x ) x , 

where Re o~(B) > 0, and relaxation systems of form 

(1.4) (») + (fa v \) =( ° 
where 

(1.5) Re a(q v (u,v*(u))) < 
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along a smooth equilibrium manifold defined by 
(1.6) q(u,v*(u)) = 0, 

and f(u) = f(u, v*(u)). Here, x, t e R 1 , u, f , f e R n , v, g, q G R r , and B G R nXn , 
where typically n, r > 1; we make the standard assumptions of strict hyperbolic- 
ity of Df, and genuine nonlinearity of the principal characteristic eigenvalue (the 
one associated with the approximate shock direction), along with further standard 
hypotheses on B and /, g, q, to be described later on: in particular, that sys- 
tems (1.3) and (1.4) be dissipative in the sense of Kawashima and Zeng [Kaw, Ze]. 
These are generically satisfied for most models of physical interest; for discussions 
of applications, we refer the reader to the general surveys [Z.3,Z.6,MZ.l,N,Yo.4]. 

A recent development in the stability analysis of smooth shock profiles has been 
the successful adaptation of Evans function/dynamical system techniques to this 
nonstandard setting, with an associated explosion of new results; see, e.g., [GZ,H.l- 
3,ZH,BSZ,BMSZ,ZS,HZ,Z.2-4,HoZ.l-2,OZ.l-2,MZ.l-2]. For the origins of these 
methods in the setting of reaction diffusion and singular perturbation problems, see 
[E.1-4,J,AGJ,PW]; see also the important analyses of [Sat,K.l-2,JGK] in the scalar 
case. In particular, it has now been established under extremely general circum- 
stances that both linearized and nonlinear stability are implied by spectral stability 
of the linearized operator about the profile, a fact that was left in doubt by earlier 
"direct" analyses [MN,Go.l-3,KMN,SX.l-2,S,GoM,LZ.l-2,FreL,L.l-3] carried out 
on a case-by-case basis. This reduces the question of stability to an ODE issue, and 
raises the hope that at least small amplitude stability can be treated in a uniform 
way across dimension and type of regularization. 

On the other hand, spectral stability has up to now been verified only by energy 
estimates resembling those of the earliest direct analyses of [MN,Go.l-2,KMN]; see 
[Z.l,HuZ,Liu,Hu.l,Z.4] for examples in various contexts. And, these have proven 
difficult to generalize: in particular, one-dimensional stability of general small am- 
plitude relaxation profiles, and multi-dimensional stability of small amplitude vis- 
cous shock profiles remain outstanding open problems in the theory. In this, and a 
companion paper [PZ] treating the multidimensional case, we remove this trouble- 
some gap in the theory, showing that the same Evans function tools that were used 
to effect the reduction to spectral stability may be used as the basis of a general 
procedure to verify spectral stability in the zero-shock-amplitude limit, applicable 
in particular in the two open cases mentioned just above. This both unifies and 
completes existing stability theory, at least as regards the small amplitude case; 
stability of large amplitude shock profiles remains as the outstanding open question 
in this area. 

The basic idea behind our approach is a natural one suggested by Gardner and 
Jones [GJ.3] in the context of the one-dimensional strictly parabolic case. In the 
weighted energy method of Goodman [Go. 1-2], the approach on which most sub- 
sequent extensions have been based, the strict transversality of characteristic fields 
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other than the principal one is used to "project out" behavior in transverse modes, 
reducing the problem to an approximate scalar conservation law in the principal 
characteristic field. Gardner and Jones pointed out that transverse modes corre- 
spond in the associated eigenvalue ODE roughly to "fast" and "super-slow" modes, 
which on a large portion of the relevant spectral domain may be projected out 
using dynamical system techniques to leave an approximately scalar "slow" man- 
ifold flow as in the argument of Goodman, and proposed this observation as the 
basis of a dynamical systems argument alternative to the one-dimensional stability 
argument of Goodman. However, Gardner and Jones did not examine the crucial 
small- frequency regime where super-slow and slow modes intermingle, and therefore 
eliminated only the possibility of eigenvalues with real part greater than a small, but 
fixed constant 9 > 0; by contrast, a standard Garding-type energy estimate gives 
the much stronger result Re A < Ce 2 , Im A < Ce, where e := \u+ — it_| denotes 
shock amplitude. Thus, at a technical level, substantial issues still remain; indeed, 
virtually all of our effort will be directed to the understanding of the "mixing" 
regime complementary to that considered in [GJ.3]. 

In the remainder of the introduction, we give an overview of the analysis and 
main results, deferring technical aspects to the main body of the paper. 

Spectral stability. Let us first make precise the notion of spectral stability 
that we consider here. Without loss of generality taking s = in (1.1) (by shifting 
to an appropriate traveling frame), we have that u(x) is a stationary solution of the 
nonlinear evolution system 

(1.7) u t = F{u) 

described by its associated regularized conservation law, of form (1.3) or (1.4). 
Linearizing this system about the stationary solution u, we obtain a linearized 
evolution system 

(1.8) v t = Lv, 

where L := dJ r (u) is the formal derivative of operator T about u. 

Definition 1.1. Following [ZH,Z.3], we define strong spectral stability (condition 
(Dl') of [Z.3]) as 

(1.9) Op(L) C {A : Re A < 0} U {0}. 

where a p (M) denotes point spectrum of a linear operator M: equivalently, there 
exist no L 2 solutions of the eigenvalue equation (L — X)w = for Re A > and 
A ^ 0. 

Remark 1.2. It has been shown in a variety of contexts, in particular, all con- 
texts considered here, that strong spectral stability together with the "hyperbolic" 
conditions of inviscid stability of the corresponding ideal shock (u-,u+) of (1.2) 
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(dynamical, or "outer" stability) and transversality of the traveling wave connec- 
tion (structural, or "inner" stability) implies linearized and nonlinear stability; see 
[ZH,Z.2,Z.3,MZ.l-2]. In the small amplitude case, under reasonable assumptions, 
inviscid stability holds always; see [M.l-3] in the one-dimensional case, [Me.l- 
5,FM,FZ,Z.5-6] in the multidimensional case. Likewise, under reasonable assump- 
tions, connections are always transverse; see discussion of existence theory given 
below. Thus, strong spectral stability is sufficient to imply linearized and nonlinear 
stability. 

The canonical model. In the one-dimensional genuinely nonlinear context, it 
is well known that Burgers equation, 



approximately describes small-amplitude viscous behavior in the principal charac- 
teristic mode. In particular, the family of exact solutions 



give an asymptotic description of the structure of weak viscous shock profiles in 
the principal direction, in the limit as amplitude e goes to zero. Note that (1.10) 
is invariant under the parabolic scaling x — > ex, t — > e 2 t, u — > u/e, whence we may 
immediately conclude that any unstable spectrum of the linearized operator about 
the profile u e must lie within a ball of radius Ce 2 . (By standard considerations, the 
unstable spectrum of a second order parabolic operator is at least bounded) . 

Burgers profiles are strongly spectrally stable, as may be readily established 
by a number of different techniques: for example, by L 1 contraction [Kr], Sturm- 
Liouville considerations [Sat, He, Z. 5], or direct energy estimates as in [MN,Go.l- 
2,Z.3 Appendix A. 6)]. We record this fact as: 

Proposition 1.3. Shock profiles (1.11) (of any amplitude e, by scale-invariance) 
are strongly spectrally stable as solutions of (1.10). 

The reduced profile equations. Existence/structure of small-amplitude pro- 
files may be determined by center manifold reduction, as pointed out by Majda 
and Pego [MP] in their pioneering treatment of the general strictly parabolic case; 
extensions to general partially parabolic, or "real" viscosity, and relaxation systems 
have been carried out in [P], [Yo.Z,FZe,MZ.l], respectively. The result in all cases, 
is that the center manifold associated with the weak profile problem (appropriately 
chosen) is foliated by one-dimensional fibers, on which the reduced flow, in the 
genuinely nonlinear case, is approximately that of the profile equation for Burgers 
equation (1.10). 

More precisely, the center manifold lies tangent to (i.e., within angle (9(e) of) 
the principal eigendirection r p of Df 1 , and the reduced flow, parametrized by 



(1.10) 



u t + (u 2 /2) x = u xx , 



(1.11) 



u e (x) := — etanh(ea;/2) 
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rj := l p (u-) ■ (u — (u- + u + )/2), l p the corresponding left eigendirection, obeys 
the approximate Burgers profile equation 

Pv' = ^(V 2 ~ e 2 ) + 0(\ V , 6| 3 ), 

or, rescaling by 

(1.12) r)^r}je, x^Aex//3: 



(1.13) V '=^(r, 2 -l) + eO(r,,e), 
0(-, ■) G C 1 , as compared to the exact profile equation 

(1.14) ff = \{f - 1) 

for (1.10) (1.11) with e = 1; here and below, " / " denotes d/dx. The effective 
diffusion coefficient (3 is given by (3 := l p Br p (u-) in the viscous case, and (3 := 
l p B*r p (ii-) in the relaxation case, where 



(1.15) 



B*(u) := -f v q v 1 {g u ~ v*J u ) 

= -fvq^ifju - 9vqy 1 q u + q^Quifu - fvqy 1 q u ), 



g := g(u, v*(u)), denotes the effective "Chapman-Enskog" viscosity associated 
with the relaxation system, and A := /*D 2 / 1 (r p , r p )(u-) ^ denotes the co- 
efficient of genuine nonlinearity for the hyperbolic system (1.2); for details, see 
[MP,Z.l,P,Yo.Z,FZe,MZ.l]. 

From this reduction to normal form we obtain in standard fashion: 

Proposition 1.4. Under the rescaling (1.12), there hold 

(1.16) \v±~V±\<Ce 
and 

(1.17) \(v-V±)-(v-V±)\<Cee- e ^ 

for x ^ for any fixed < 9 < 1, and some C = C(6) > 0, where rj± and fj± := ±1 
denote the rest points of (1.13) and (1.14), respectively. 

Proof. Assertion (1.16) follows by the Implicit Function Theorem from normal 
hyperbolicity of fj±. Assertion (1.17) follows from standard stable/ustable manifolds 
estimates. □ 
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Corollaryl.5. Under the rescaling (1.12), the principle eigenvalue a p := a p (u) 
converges to fj at the same rates (1.16)-(1.17) as does r\. 

Proof. This follows using tangency of the profile to direction r p together with 
the fact [Sm] that Va p ■ r p = A. □ 

The reduced eigenvalue equations. Following [GJ.3], we carry out here a 
reduction of the generalized eigenvalue equations similar in spirit to that carried 
out by Majda&Pego and others for the profile equations. However, whereas those 
anlayses involved center manifold reduction of an autonomous nonlinear system, we 
shall work directly with the linear, nonautonomous system given by the eigenvalue 
equation (L — X)w = associated with (1.8), written as an appropriate first-order 
system 



using the tracking lemma of [GZ,ZH] (specifically, a refined version introduced in 
[MZ.l]) to effect the reduction to a slow manifold; see Section 2. Here, (1.18) 
is dimension 2n in the strictly parabolic viscous case, n + RankS in the general 
viscous case, and n + r in the relaxation above, '"" denotes d/dx. For 

the efficient coordinatization of the first-order system (1.18) in the relaxation and 
general viscous case (important for practical computation), we follow the scheme of 
[Z.3]; see sections 3-4 for details. 

Of course the tracking lemma is itself an analytic formulation of classical invariant 
manifold techniques; see the original formulation in [GZ] in terms of projectivized 
flow/invariant sets of autonomous systems. The calculation of the reduced flow 
on the slow manifold can equally well be carried out by classical center manifold 
reduction; see the concurrent treatment of Plaza [PI] in the one-dimensional case; 
see also the independent treatment of [FSz] (described further in the note following 
the introduction). Indeed, this may be preferable in the more complicated situations 
arising with various types of degeneracies in the underlying equations. However, 
there is some advantage in taking account of the linear nature of the underlying 
problem, in that the calculation of reduced equations simplifies considerably. 

First reduction. A first application of the tracking lemma projects out "fast" 
transverse modes: n — 1 in the strictly parabolic viscous case, RankS — 1 in the 
general viscous case, and r — 1 in the relaxation case, in each case reducing to an (n+ 
1) -dimensional "slow manifold," of which n— 1 correspond to "super slow" transverse 
modes pj, and the remaining 2 to "medium slow" Burgers modes which are roughly 
linearized versions of n, rj' =: z in the limiting equations above. Rescaling by 
x — > Aex/P, A — > /3A/A 2 e 2 , and (balancing by) rj' — > firf/Ae, we obtain a normal 
form for the reduced flow that has bounded (0(1)) coefficients for bounded A and, 
modulo error terms &j(x, X)W (W := (p, 77, z)), becomes a block triangular system 
in p = (p + , p_) and (77, z), of form 



(1.18) 



W = A(x,X)W, 



(1.19) 



p'± = M±(x,\)p±, 
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and 

(1-20) (jty = M (x,X)^+N(x,X)p, 
where: 

(1.21) $ p± (±oo, A) = 0, <& (r ^)(±oo, A) = e£±(X, e) 

(1.22) |$(x,A)-$(±cx),A)| <C(l + |A|)ee- e|a; l for x ^ 0, 
and 

(1.23) \N\<Ce- e W; 

(1.24) |M±| = £>(e|A|), M± ^ e#(|Re A| + e 2 |Im A| 2 ), 

6* > 0; and M in the viscous case corresponds exactly to the coefficient 



(1.25) A := 



1 

A fj(x) 



of the eigenvalue equation for the standard Burgers equation (1.10) (1.11) , e = 1, 
written as a first-order system. Here, £±(-) is a smooth function satisfying the 
growth bounds 

ri2fi N /(9(l + |A|e) 0(l + |A|e)\ 

{1 - 2b) 6± - v O(l + |A|) 0(1 + |A|);- 

Remark 1.6. The asymmetric bounds on S± are a result of rescaling. The stronger 
bound in the upper righthand corner is actually important for our later analysis, 
specifically, in the large-frequency regime where we approximately undo the rescal- 
ing in the course of our argument. 

Further reductions/normal forms. In the rescaled coordinates, there are three 
distinct frequency regimes: 

I. |A| < C; II. C < |A| < Ce" 1 ; and III. Ce" 1 < |A| < C"^" 2 , 

C > sufficiently large, corresponding in the original (unrescaled) coordinates to 
|A| < Ce 2 , Ce 2 < |A| < Ce, and Ce < |A| < C -1 , respectively. (By standard 
considerations, the spectrum of L is restricted to small frequencies |A| < C _1 ; see 
Proposition 2.9, below.) Each regime requires a slightly different treatment, and is 
associated with a different normal form. 
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Region III, the region considered by Gardner and Jones [GJ], we will denote as 
the "parabolic" regime, as this is the regime on which behavior is dominated by 
dissipativity of the underlying system (which yields "parabolic" behavior for small 
frequencies, as evidenced by the structure of the reduced system (1.19)-(1.20)), 
independent of any other structure of the traveling wave profile whatsoever. This 
regime may be treated by the standard "high-frequency" techniques of the usual 
Evans function theory (see [AGJ,GZ,ZH,MZ.l]), which are roughly equivalent to 
sectorial, or Garding-type, energy estimates, to preclude altogether the possibility 
of eigenvalues. 

Specifically, (r),z) may be resolved into decaying/growing modes rj± that have 
a uniform spectral gap both from each other and from p + and p_, the latter of 
which also have a uniform gap from one another, to an order sufficient that a 
second application of the tracking lemma, to the once-reduced equations (1.19)- 
(1.20) yields further reductions of the flow onto four transverse invariant manifolds 
tangent to each of the p±, rj± directions. In particular, we find that the manifold 
of solutions decaying at +oo is transverse to that of solutions decaying at — oo, 
yielding the result. In other words, the normal form for this regime is trivial, with 
growing and decaying modes decoupled. 

Region II we denote as the "reduced parabolic" regime, since this is the regime for 
which similar tracking/energy estimates prohibit eigenvalues for the Burgers equa- 
tion. Here also, we find that eigenvalues cannot exist, independent of the structure 
of the profile. However, the analysis is much more delicate; indeed, we regard this 
as the trickiest case in our argument. Here, we can again separate off decay/growth 
modes rj± (parabolic behavior of the Burgers part) as in the previous case, but p± 
are not sufficiently separated to apply the tracking lemma to these coordinates. In 
this case, therefore, our ultimate reduction is to the (n — l)-dimensional superslow 
equations (1.19), but with error terms now involving p alone and not (rj,z). Note 
that the coefficients M± are bounded on regime II, so that this is indeed a valid 
normal form. 

Region I we denote as the "gap regime," in reference to the fact that on this 
regime we rely solely on the gap lemma of [GZ,KS,ZH,Z.3] (specifically a refined 
version given in [MeZ]) in our analysis of the slow flow, rather than using the 
tracking lemma to carry out a further reduction as in Regions II and III. The "gap" 
in "gap lemma" refers in fact to absence of spectral gap, as is the case in this 
regime; the lemma asserts that uniform exponential convergence of coefficients as 
x — > ±oo and continuous extensionof the stable/unstable eigenspaces of the limiting 
coefficient matrices A±(A) at x — > ±oo can substitute in this situation for uniform 
spectral gap, to yield estimates valid on x ^ 0. In this regime, which is the crucial 
one determining behavior, the normal form is the entire (n + l)-dimensional system 
(1.19)-(1.20). 

The gap lemma: convergence to limiting flows. We complete our reduc- 
tion in each case by an application of the gap lemma (specifically, a refinement 
given in [MeZ]), showing that the flow of the normal forms we have obtained con- 
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verges to that of the "limiting" normal forms obtained by dropping error terms 
OW, etc.: more precisely, the flows associated with the stable manifold at +00 
and the unstable manifold at —00, or in some cases their analytic continuation into 
the essential spectrum. The basic approach is to: (i) conjugate approximate and 
limiting equations on the half-lines x ^ to their asymptotic systems at ±00 by a 
linear change of variables uniformly close to the identity (an application of the gap 
lemma, made possible by the uniform exponential decay estimate (1.21)), then (ii) 
check by direct linear algebraic computation that the stable and unstable subspaces 
of the respective asymptotic coefficient matrices are uniformly close. In the present 
applications, (ii) is straightforward, since the asymptotic matrices of approximate 
and limiting flows agree for the p equation, while for the (decoupled) (77, z) equa- 
tion, the stable and unstable subspaces for the limiting equations are, under our 
strategically chosen rescaling, uniformly spectrally separated on the region I where 
it comes into play. 

Convergence of the stable flow at +00 and unstable flow at —00 then implies 
uniform convergence as e — > of suitably chosen Evans functions associated with 
the approximate and limiting normal forms. The Evans function is a Wronskian 
measuring solid angle between the stable and unstable subspaces of a given eigen- 
value ODE written as a first-order system, constructed so as to be analytic in the 
frequency A on a "region of consistent splitting" (defined (Al), Section 2). In each 
of the contexts considered here, the region of consistent splitting includes the en- 
tire set of frequencies {Re A > 0} \ {0} of interest; see, e.g., [Z.3]. Moreover, both 
the Evans function and its component subspaces extend continuously along rays in 
frequency space, in both the one- and multidimensional case [ZS,Z.3] to the closure 
at the origin, making possible uniform estimates of transversality; indeed, in the 
one-dimensional case considered here, the Evans function may in fact be extended 
analytically through the origin, without restricting to rays [GZ]. (Note: in the mul- 
tidimensional case, the Evans function may be holomorphically extended along rays 
[ZS,Z.3]; however, in certain "glancing" directions, there arise branch singularities 
at the origin.) On the region of consistent splitting, the zeroes of the Evans function 
correspond both in location and multiplicity with the eigenvalues of the associated 
linear operator; for details, see [GJ.1-2,ZH]. 

The multiplicity of the generalized eigenvalue at the origin A = is directly calcu- 
lable and related to hyperbolic stability, as shown in the one- and multi-dimensional 
case in [GZ,BSZ] and [ZS], respectively. In the present, small-amplitude case, it is 
readily seen to be the same as the multiplicity for the limiting normal form equa- 
tions, namely one. In the one-dimensional, analytically extendable case, we may 
thus conclude immediately from uniform convergence to the limiting Evans func- 
tion (by degree, i.e., winding number, considerations) that strong spectral stability 
of the original system is equivalent to strong spectral stability for the limiting normal 
forms. Alternatively, we may remove the zero of the Evans function at the origin 
by working instead with the "integrated" eigenvalue equations following [MN,Go.l] 
or with "flux variables" following [Go. 2], to recover the same conclusion using only 
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continuity at the origin of D. 

In this paper, we follow the latter strategy exclusively. In the multidimensional 
analysis of [PZ], we use a modification of the same strategy, specifically, an inter- 
esting variation intermediate to flux and integrated variable methods, neither of 
which themselves directly generalize to multidimensions in a useful way, to achieve 
the same result; this "balanced flux form" is quite similar to that introduced in 
the treatment of relaxation systems in Section 4 of the present paper. The former 
strategy is perhaps viable as well, using computation of winding number on Rie- 
mann surfaces as suggested in [GZ]; however, this would involve extra bookkeeping 
(specifically, tracking of branch points in approximate vs. limiting systems) that 
we prefer to avoid. 

Conclusions/summary of the main results. From this point, we may im- 
mediately deduce our main results. For, in each of the (nontrivial) limiting normal 
forms arising in regions I II, the equations for p± completely decouple, and clearly 
support no "decaying" solutions other than the trivial, zero solution (in general, 
"decaying" is defined as lying in p + direction at +00 and p- direction at —00; away 
from the essential spectrum, this is equivalent to actual decay at ±00). Thus, we 
may conclude that p = for any generalized eigensolution, leaving a decoupled 
(rj, z) equation in region I, agreeing with the generalized eigenvalue equations for 
canonical model (1.10), and the trivial flow in region II. Recall that the reduced 
flow was already trivial in region III. Likewise, v = in each of these regions, since 
fast growing and fast decaying modes completely decouple. 

Noting that the canonical model, being among the class considered, is subject 
to the same reduction in regions II III, we may conclude, finally, our main result: 
that strong spectral stability is equivalent to strong spectral stability of the canonical 
model, at least modulo the degenerate case that the canonical model yields nontrivial 
nonstable eigenvalues lying precisely on the imaginary axis; see Proposition 3.2 and 
Remark 3.3 below. As a corollary, recalling the result of Proposition 1.3, we obtain 
the stated results of strong spectral stability of small-amplitude profiles, across the 
general class of models considered in (1.3) and (1.4). More precisely, in this paper, 
we carry out the details of the one-dimensional case for general strictly parabolic 
viscosities, and for general relaxation systems. The case of general real viscosity 
may be carried out similarly as the relaxation case, following the dual framework 
set out in Appendix A. 2 of [Z.3]; we omit these calculations, as satisfactory one- 
dimensional results already exist [HuZ]. 

Plan of the paper. In Section 2, we recall and slightly extend the basic Evans 
function tools we will need. In Section 3, we carry out the one-dimensional, strictly 
parabolic viscous case. For clarity, we first carry out in its entirety the simpler 
identity viscosity case, afterward indicating the necessary adjustments in the general 
case. In Section 4, we carry out the general relaxation case, for clarity first treating 
the 2x2 case by reduction to the scalar viscous case; we note that this gives an 
independent proof of stability in this case, different from that of Liu [L.2]. Finally, in 
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Section 5, we briefly indicate the extension of our analysis to the multidimensional 
viscous case; details will be given in [PZ]. 

Note: Shortly before the completion of this analysis (precisely: after the comple- 
tion of the one-dimensional viscous case and before the completion of the relaxation 
case), we have learned of closely related work of H. Freistiihler and P. Szmolyan 
[FreS] in which they establish a reduction method similar (indeed, perhaps equiva- 
lent) to ours, but proceeding entirely from the point of view of geometric singular 
perturbation theory. In the cited work, the first in a planned series of three, they 
carry out a complete analysis of the one-dimensional identity viscosity case, an- 
nouncing the intention to treat one-dimensional relaxation/real viscosity and mul- 
tidimensional stability in papers two and three, respectively. 

Section 2. Evans function framework. 

We begin our analysis by recalling, and in some cases slightly extending, the 
Evans function tools we will need. For later reference, we carry out the analysis in 
generality sufficient for the multidimensional case as well. 

2.1. The gap lemma and convergence of approximate flows. Consider a 
family of first-order systems 

(2.1) W = A e (x,l\)W 
indexed by parameter e, where (£, A) vary within the set 

(2.2) O := {(I, A) : £ G M. d ~\ Re A > 0} 

of unstable-neutrally stable frequencies, e varies within some given open set V, and 
x varies within R 1 . Equations (2.1) are to be thought of as generalized eigenvalue 
equations, with parameter A representing a Laplace transform frequency in time, 
and £ representing a Fourier transform frequency in directions of spatial symmetry; 
in the present, one-dimensional case, £ = 0. In the applications of this paper and 
of [PZ], e will be just the shock strength \u+ — U-\. 
We make the following basic assumptions. 

(AO) Coefficient A e (-,£, A), considered as a function from (£, A, e) into L°°(x) 
is analytic in (£, A) and C k in e for some k > on O x V. Moreover, A e (-,£, A) 
approach exponentially to limits A± as x — > ±oo, with uniform exponential decay 
estimates 

(2.3) |A e - A £ ± | < Cie- |x|/C2 , C, > 0, for x ^ 

on compact subsets of O x V. 

(Al) On (O\{(0, 0)}) x V, the limiting matrices A^ and Af_ are both hyperbolic 
(have no center subspace), and the dimensions of their stable (resp. unstable) 
subspaces S+ and St (resp. unstable subspaces U+ and Ut) agree. 
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(A2) At the origin (£, A) = (0,0), subpaces S±, U± have continuous (to some 
order, which may be even analytic, or holomorphic with branch point at the origin) 
limits along rays (i.e., for (f , A) = (r| , r\ ) as r G R -> 0+, (| , A ) G O \ {(0, 0)}), 
which, moreover, are C k in the parameter e, where is as in (AO) above. 

These are satisfied in all the contexts considered in both this paper and [PZ]; 
see [Z.3]. Condition (AO) is induced, at least for e bounded from zero, by the ori- 
gins of (2.1) through the linearization about smooth shock profiles of (1.3) or (1.4). 
Uniform approach (2.3) as e — > will be recovered through reduction/rescaling. 
Condition (Al) may be recognized as the standard hypothesis of "consistent split- 
ting," as introduced in [AGJ], and follows by the assumption of dissipativity of 
systems (1.3) and (1.4). Condition (A2) is an extension suitable for multidimen- 
sions of the "geometric separation" hypothesis of [GZ] and is a consequence of the 
hyperbolic structure of (1.2). (Note: separation was not needed for the arguments 
of [GZ], and in fact does not hold in multidimensions; see discussion of Appendix 
A.4, [Z.3].) 

The gap lemma. Under these circumstances, the "gap lemma" established in 
various degrees of generality in [GZ,KS,ZH,Z.3] asserts that behavior of the variable- 
coefficient equation (2.1) may be related to that of its constant-coefficient limiting 
systems on x ^ 0, while maintaining the assumed regularity in all parameters. This 
has recently been greatly improved in [MeZ] , in the form of the following conjugation 
lemma, a version that is particularly convenient in applications. 

Lemma 2.1 [MeZ]. Under assumptions (A0)-(A2), there exists locally to any given 
(£, A, e) G O x V a pair of linear transformations PI (x, £, A) and PL (x, £, A) on x > 
and x < 0, respectively, and possessing the same regularity as A e in all arguments, 
such that: 

(i) \P±\ and their inverses are uniformly bounded, with 

(2.4) \P1_ -I\< CCiC 2 e- e|x|/C2 

where < 9 < 1 is an arbitrary fixed parameter, and C > and the size of the 
neighborhood of definition depend only on 9, the modulus of the entries of A e at 
(£, A, e) G O x V, and the modulus of continuity of A e on some neighborhood of 
(f,A,e)GfixV. 

(ii) The change of coordinates W := P±Z reduces (2.1) on x > and x < 0, 
respectively, to the asymptotic constant- coefficient equations 

(2.5) Z' = M ± (l\)Z. 

Proof. As described in [MeZ] , this is a straighforward corollary of the gap lemma 
as stated in [Z.3], applied to the "lifted" matrix- valued equations for the conjugating 
matrices P±. □ 
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Definition 2.2. Let vf + , . . . ,vl + and v^ k+1 ^_, . . . , v e N _ be bases for and 
Ut, as defined in (Al), chosen with the same regularity assumed on S+, Ut (note: 
that such a choice is possible is a consequence of standard matrix perturbation 
theory [Kat]), and P± be transformations defined as above on some neighborhood 
of (£, A, e G Ox V. Then, the local Evans function for (2.1) associated with these 
choices is defined as 

(2.6) D*(iX) :=det (Plv{ + ... P-.v% + Pl^+i)- ••• P - v n-) {x=0 - 

Remark 2.3. Evidently, the choice of local Evans function is highly nonunique. 
Though we shall not need it here, it can be shown using the original formulation 
of the gap lemma in [GZ,KS] that a globally analytic choice is possible in both the 
one- and multidimensional case; see [GZ] and [ZS,Z.3], respectively. 

Combining (2.5) and (2.6), we immediately obtain the following result sufficient 
for our needs. 1 

Proposition 2.4. For (£, A) within the region of consistent splitting Q \ {(0,0)}, 
equation (2.1) admits a nontrivial solution W G L 2 (x) if and only if -D e (£, A) = 0. 

Proof. Evidently, the first K columns of the matrix on the righthand side of 
(2.6) are a basis for the stable manifold of (2.1) at x — > +oo, while the final N — K 
columns are a basis for the unstable manifold at x — ■> — oo. Thus, its determinant 
vanishes if and only if these manifolds have nontrivial intersection, and the result 
follows. □ 

Convergence of approximate flows. We now turn to the crucial question: 
Under what circumstances does the Evans function of (2.1) converge as e — > + to 
the Evans function of the limiting equations at e = 0? The following proposition 
gives a simple and sharp answer. 

Proposition 2.5. Suppose that, in an Vt-neighborhood of some (£, A): 

(i) As e — > + , the asymptotic subspaces S^_, Ut converge uniformly in angle to 
S+, U^_, with rate S(e): equivalently, for < e < eo, their spanning bases satisfy 

(2-7) \v* ± -v° ± \<r,(e). 

(ii) The coefficient matrices A e converge uniformly exponentially to their limits, 
in the sense that, for < e < eq, 

(2.8) |(A e - h±) - (A - A^)| < C lV (e)e-^ /c \ 



1 Recall that Gardner and Jones [GJ.1-2] have established the stronger result of correspondence 
up to multiplicity between eigenvalues and zeroes of D. 
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Then, on some (possibly smaller) Q-neighborhood of (£, A), the local Evans func- 
tion D e defined as in Definition 2.2 converges uniformly as e — > + to D° , with 

(2.9) \D e -D°\<CC 1 C 2 v(€). 

(More precisely, the columns of the defining matrix on the righthand side of (2.6) 
converge uniformly, with the same rate CC\Cir\{e).) 

Proof. Clearly, it is sufficient to prove the stronger parenthetical assertion of 
convergence of individual columns, and for this we may restrict without loss of 
generality to the + columns, and the right half-line x > 0. Applying the conjugating 
transformation W — > {P^_)~ X W for the e = equations, we may reduce to the case 
that A is constant, and P° = 0. In this case, (2.8) becomes just 

\(A*-Al)\<C lV (6)e-W/ c >, 

and we obtain directly from the conjugation lemma, Lemma 2.1, the estimate 

\P+ -P+\ = \P+-I\< CCiC 2 ?7(e)e- e|x|/C2 , 

and in particular 

(2.10) |P;-P°|| X=0 < CdCs^e). 

The result now follows, from (2.6), (2.7), and (2.10). □ 

Remark. The same argument shows that the proposition is sharp, since easy 
examples show that the conjugation lemma itself is sharp. 

2.2. The tracking lemma and reduction of the flow. Next, consider com- 
plementary situation of a family of equations of form (2.1) on an e-neighborhood for 
which the coefficient A e does not exhibit uniform exponential decay to its asymp- 
totic limits, but instead is slowly varying. This occurs quite generally for rescaled 
eigenvalue equations arising in the study of the large- frequency regime; see, e.g., 
[GZ,ZH,MZ.1,Z.3]. In the present context, it arises for the unrescaled equations in 
the small-shock strength limit e — > 0. 

In this situation, it frequently occurs that not only A e but also certain of its 
invariant (group) eigenspaces are slowly varying with x, i.e., there exist matrices 

(2.11) L*=^(x), W = {R\ Rl){x) 
for which L e R e (x) = I and \LR'\ = \L'R\ is small relative to 

(2.12) M £ :=L e A^(x)= ^(x), 
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where '"" as usual denotes d/dx. In this case, making the change of coordinates 
W e = R e Z, we may reduce (2.1) to the approximately block-diagonal equation 

(2.13) Z e ' = M £ Z e + 5 e G e Z e , 

where M e is as in (2.12), <d e (x) are uniformly bounded matrices, and 5 £ (x) < 5(e) is 
a (relatively) small scalar. A sometimes crucial improvement may be obtained by 
arranging that error 6 vanish on diagonal blocks, i.e., that LjR'- = 0; see [Go.l- 
2,MZ.l]. Here, we shall make use of this observation only in the principal, 2x2 
"Burgers" block. 

Let us assume that such a procedure has been successfully carried out, and, 
moreover, that there exists a uniform spectral gap, in the strong sense that 

(2.14) Re (Mi - M 2 ) > 77(e) > Ofor all x. 

In the "standard" case that M e are uniformly bounded and 77(e) may be taken inde- 
pendent of e, that (2.14) may be arranged by a suitable coordinate transformation 
provided that there holds the weaker condition 

(2.15) minRe a(M 1 ) - max Re a{M 2 ) > fj > Ofor all x, 

for some 77 > 77. Then, there holds the following reduction lemma, a refinement 
established in [MZ.l] of the "tracking lemma" given in varying degrees of generality 
in [GZ,ZH,Z.3]. The new feature of the reduction lemma is that it asserts the 
existence of smooth invariant manifolds in the vicinity of the first and second block 
coordinates, whereas the tracking lemma(s) assert only the existence of forward and 
backward attracting cones in the same vicinity and do not give a reduction in the 
usual dynamical systems sense. 

Proposition 2.6[MZ.l]. Consider a system (2.13) under the gap assumption 
(2.14), with O e uniformly bounded for e sufficiently small. If, for < e < 6q, 
the ratio 5(e) /77(e) of formal error vs. gap is sufficiently small relative to the bounds 
on Q, in particular if S(e)/rj(e) — > 0, then, for all < e < eo, there exist (unique) 
linear transformations 3>\(x, £, A) and $|(x,^, A), possessing the same regularity 
with respect to the various parameters e, x, A as do coefficients MT and 6(e) e , 
for which the graphs {(Z\, $2-^1)} an d {(^1(^2); are invariant under the flow 
of (2.13), and satisfying 

\<5>\\, \ < C8(e) /77(e) for all x. 

Proof. This may be established by a contraction mapping argument, carried 
out for the projectivized "lifted" equations governing the flow of exterior forms; see 
Appendix A. 3. 2. 2 (ii) of [MZ.l]. See, e.g., [Sat] for a corresponding argument in 
the case that Mi, M 2 are scalar; the lifting to exterior forms essentially reduces the 
problem to this case. □ 
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Prom Proposition 2.6, we obtain reduced flows 

(2.16) Z{' = M{Z{ + 8 e e\ 2 Q\Z{ = M\Z\ + 0(5 e {x))Z{ 
and 

(2.17) Z|' = M|Z| + ^eii^ZJ = MJZJ + 0(5 £ (x))Z| 

on the two invariant manifolds described. Should we arrange that ©n = or ©22 = 
0, then the error terms would become 0(S e (x)5(e) /rj(e)) Z{ or (D(5 e (x)5(c) /rj(c)) Z%, 
respectively; that is, enforcing vanishing of on a diagonal block reduces the error 
term by factor 5/rj. 

Remark 2.7. As pointed out in [MZ.l], the gap condition (2.14) may be allowed 
to fail by order a on an interval of length order 1/a, for any fixed a. For the 
applications of this paper and of [PZ], as for those of [MZ.l], this means that the 
gap condition need only be checked at x = ±00. 

Remark 2.8. As described in [MZ.l], provided that Mi and M2 are (uniformly) 
bounded and spectrally separated, one can in fact repeat the procedure used to 
construct (2.13) to obtain an approximate block-diagonalization to arbitrarily high 
order: i.e., an asymptotic expansion in powers of 5(e). In general, e.g., in the 
large-amplitude situation considered in [MZ.l], these higher order terms are not 
explicitly computable. However, we point out that, in the present small-amplitude 
situation, they may in principle be determined exactly to any desired order, using 
the error expansion for the traveling wave profile to approximate L'; indeed, this 
can be recognized as just a variation of the usual center manifold calculations used 
to approximate center manifold flow. In the calculations of this paper and of [PZ], 
we shall not need to perform any such higher-order corrections: a single iteration 
will be enough. 

2.3. Reduction to small frequency. As a first application of Proposition 
2.6, we obtain the following preliminary result analogous to that of [GJ], reducing 
the small-amplitude stability problem to the small frequency, or "diffusive," regime 
|A| « 1. 

Proposition 2.9. In each of the contexts considered in this paper ( one- dimensional 
viscous and relaxation systems, with hypotheses given below), the nonstable eigenval- 
ues Re A > of linearized operator L of (1.8) are restricted to |A| < r(e), r(e) — > 
as e — > 0, where e := \u+ — U-\ denotes shock strength. 

Proof. Equivalently, we show, for arbitrary r > 0, that |A| < r for e sufficiently 
small. For |A| > r and A e as in (2.1) denoting the coefficient of the generalized 
eigenvalue problem (1.8) written as a first-order system, there is a uniform spectral 
gap between stable and unstable subpaces of A e (x), independent of < e < eo and 
x; as discussed in [Z.3], this is a straightforward consequence of the assumed dissi- 
pativity of the systems (in the sense of Kawashima-Zeng [Kaw,Ze]). By standard 



RAMON PLAZA AND KEVIN ZUMBRUN 



17 



matrix perturbation theory [Kat] plus boundedness (compactness) of A e , we may 
therefore deduce the existence of well-conditioned bases (2.11) for these subspaces 
varying smoothly with A e , and therefore satisfying 

\LR'\ = 0(\(^y\) = 0(e 2 ) 

by Proposition 1.4. Moreover, by a proper choice of basis, we may arrange (2.14) 
with r\ > uniformly bounded, independent of < e < eo and x; for details, see 
Appendix A4, [Z.3]. 

Applying Proposition 2.6, we obtain a pair of reduced systems (2.16) and (2.17) 
for which the coefficent matrices (with 0(e 2 ) error term taken into account) are, 
respectively, uniformly positive and negative definite, and solutions therefore uni- 
formly exponentially grow and decay, provided that e is sufficiently small. It follows 
that the equations have no solutions bounded at plus and minus spatial infinity save 
for the trivial solution W = 0. 

For |A| > R, R > sufficiently large, on the other hand, the nonexistence of 
eigenvalues may be established independent of the size of e by a rescaling followed 
by a similar reduction argument; see [ZH,MZ.l-2]. for a multidimensional version, 
see Lemma 4.38, [Z.3]. □ 

Remark 2.10. The same argument applies word for word in the multidimen- 
sional case, to yield |(£, A)| << 1, where £ denotes the Fourier transform frequency 
in directions parallel to the shock front, and A as here denotes the Laplace transform 
frequency with respect to time. 

Section 3. One-dimensional stability for parabolic systems. 

With these preparations, we are ready to carry out the stability analysis. We be- 
gin with a treatment of small-amplitude stationary profiles u £ (x) of one-dimensional 
strictly parabolic viscous conservation laws 

(3.1) u t + f(u) x = {B(u)u x ) x 

in a neighborhood U of a particular state u$. We make the assumptions: 

(HO) f,BeC 2 (regularity). 

(HI) Re a(B(u)) > (strict parabolicity) . 

(H2) a(Df) real, simple (strict hyperbolicity). 

(H3) Re a(-iDfC - B^ 2 ) < -#|£| 2 , f € R, for some 6 > (linearized stability 
of constant states). 

(H4) a p = is a simple eigenvalue of Df(u ) with left and right eigenvectors l p 
and r p , and l p D 2 f(r p ,r p ) ^ (genuine nonlinearity of the principal characteristic 
field). 

Here, e > denotes shock strength \u e + — u € _\, and profiles u e (-) are assumed to 
converge as e — > to uq. Under these assumptions, the center-manifold argument 
of Majda&Pego [MP] verifies the assertion of Proposition 1.4, yielding convergence 
after rescaling of to the standard Burgers profile (1.11). In the rest of this section, 
we shall establish our first main result: 
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Theorem 3.1. Under assumptions (H0)-(H4), profiles u e are strongly spectrally 
stable (and therefore linearly and nonlinearly orbitally stable [ZH,Z.2]) for e suffi- 
ciently small. 

More precisely, we shall establish the following proposition, from which, together 
with Proposition 1.3, Theorem 3.1 follows as an immediate corollary. 

Proposition 3.2. Under assumptions (H0)-(H4 ), profiles u e are strongly spectrally 
stable (and therefore linearly and nonlinearly orbitally stable [ZH,Z.2]) for e suffi- 
ciently small, if the standard Burgers profile (1.11) is strongly spectrally stable for 
e = 1 and only if the Burgers profile is not strongly spectrally unstable for e = 1 in 
the sense that there exists an eigenvalue Re A > 0. 

Remark 3.3. Of course, the "only if" part of Proposition 3.2 follows vacuously 
in the case at hand; the point is that we obtain this information by our method of 
proof, independent of any knowledge of behavior of the the limiting system. This 
distinction may be important in future applications to more general situations. 

Remark 3.4. Theorem 3.1 recovers the spectral stability result obtained by 
energy methods in [Go. 1-2], and slightly extends it from the class of A, B such 
that LBR > for some diagonalizing bases L, R = L~ x of A to the class of all 
"stable" pairs A, B in the sense of Majda&Pego [MP]. (In the case that A, B 
are simultaneously symmetrizable, these two classes are equivalent; see [MP].) It 
is straightforward with some additional bookkeeping (see [Z.3,HuZ]) to extend our 
results to the case, to which the nonlinear analyses of [ZH,Z.2] also apply, that (H2) 
is relaxed to: 

(H2') cr(Df) real, semisimple. 

3.1. The case B = I. To illustrate the method, we first carry out the proof 
of Proposition 3.2 and Theorem 3.1 in the setting of identity viscosity B = I, for 
which the associated linear algebra is particularly simple. 

Linearizing about u e , we obtain the family of eigenvalue equations 

(3.2) w" = (A e w)' + Xw, 
where 

(3.3) A^(x) := Df(u*(x)). 

Following Goodman [Go.l], we consider instead the more favorable "integrated 
equations" 

(3.4) w" = A e w' + Xw 

satisfied by the anti-derivative w(x) := f*^ w of w. For A G {Re A > 0} \ {0}, it is 
easily seen that (3.2) possesses I? solutions if and only if does (3.4) [ZH]: indeed, 
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the L 2 solutions of (3.2) consist exactly of the derivatives of those of (3.4). Unlike 
(3.2), which has a single zero eigenfunction w = (u £ )' , (3.4) possesses no zero (L 2 , 
i.e., decaying) eigenfunctions. For, the only candidate, w = u e — u e _, does not decay 
as x — > +oo, but converges to ul — u e _. 

Remark 3.5. It is worth noting that equations (3.4) may be obtained by another 
route. As pointed out by Goodman [Go. 2], the flux variable w := w' — Aw satisfies 
the same equation. 

Writing (3.4) as the family of first-order systems 

(3.5) W = A e (x,X)W, 
where 

(3.6) W :=(:,), A.:=g [. ) , 
we are ready to begin our analysis. 

Lemma 3.6. The eigenvalues fx and associated left and right eigenvectors L, R of 
A e (x) in (3.6) may be denoted as 

aj ± Jo 2 - +4A 

(3-7) 4(x) = V_ ; 

(3-8) Lf(x) = (A/ J /(A + ^ ±2 ),^/ J /(A + ^ 2 )), Rf(x) = (^;.) , 

where aj(x) denote the eigenvalues of A e (x) and l{x), r(x) associated left and right 
eigenvectors, normalized so that Ijrt = S J k (here, we have suppressed the e depen- 
dence for the sake of readability) . 

By (H2) and (H4), we may order the (real) eigenvalues of A e as 



(3.9) ai,...,a p _i < -6 < 0, 

(3.10) a p ~ ef](Aex/P = 0(e), 
and 

(3.11) < 9 < a p+ i,...,a n , 



where A and [3 denote the genuine nonlinearity and effective diffusion coefficients 
described in the introduction, and fj the standard Burgers profile (1.10). Like the 



20 



STABILITY OF SMALL AMPLITUDE SHOCK PROFILES 



principal eigenvalue a p , the transverse eigenvalues aj, j 7^ p vary within an 0(e) 
neighborhood of the corresponding eigenvalues of Df{uo). 

Evidently, for the indices j 7^ p such that aj is bounded from zero, the eigenvalues 
are jointly analytic in aj and thus A e and A within neighborhoods of Df(uo), 
0; a brief calculation shows that and are bounded, hence analytic in these 
parameters as well. The two-dimensional group eigenspace spanned by Lp and 
Rp likewise varies analytically as subspaces, as A e and A are varied; however, the 
individual eigenvectors that make it up do not (indeed, blow up for a p = and 
A — > 0). We therefore do not attempt to resolve these subspaces, but just choose 
convenient analytically varying spanning bases 



(3.12) 
and 

(3.13) 



(l p 0),L Z :=(0 l p ), 



R v :— 







R* := 







Summing up, we have: 
Lemma 3.7. In the transverse fields j 7^ p, there hold expansions 



(3.14) 

for j < p and 
(3.15) 



fj+ = -X/aj + X 2 /a 3 J + ... 



aj + 



fif = aj + ..., 



H = -X/aj + X 2 /a 3 j + ... 



for j > p, analytic in X and smooth (indeed, analytic) in aj, or, through aj, in A e , 
on sufficiently small neighborhoods of X = and aj = aj(uo) or A e = Df(ua). Like- 
wise, the eigenvectors , R^ , j ^ p possess the same regularity, as do bases L v , 
L z and R v , R z (defined in (3. 12) -(3. 13)) for their complementary two-dimensional 
invariant subspaces. 

Proof. Taylor expansion/direct calculation. □ 

Approximate block-diagonalization. Following the procedure outlined in 
Section 2.2, define now 

( v ~\ 

P- 

V 

z 

P+ 



(3.16) 



Z = 



:= L e W, 
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(3.17) 



where 



(3.18) 



(3.19) 



(3.20) 



/ L u _ 




L P - 








\L U+ 




) 



L v _ := 



(Rv_ , , Ro, Ru + j • 



-R^_ : — (^Ri , • • • , 1 ' ) ' 



-R p _ :— (^Rp + i, . . . , R^, ^ j 



(3.21) 



and 



(3.22) 



( Lr 

V L~ ) 



Rp+ '■— (R-i 7 • • • 7 R p -n ) ? 



R v+ ■— (^R p+ i, ■ ■ • , R n ^ , 



L e M e = i", where all coefficients depend on x through A e (u e (x)). 

This is just such an approximately block-diagonalizing transformation as de- 
scribed in Section 2.2: in particular, because L e , IR e depend smoothly on A e , we 
have that 



(3.23) 



6 £ := Ir/r^W = 0(1), 



where rj := l p (u-)-(u e — (ii-+u+)/2) as in the introduction. Thus, in Z coordinates, 
(3.5) reduces to the form (2.13), with 



(3.24) 



8'(x) := \rj'\ = 0(e 2 e- e ^), 6 > 0, 
5(e) = 0(e 2 ), 
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(3.25) |e e |=C(l), 
and 

(3.26) M e := diag(M u _,M p -,M ,M p+ ,M v+ ), 
where, for Re A > 0: 

(3.27) M v _ : = diag (//+,..., ^-i) < ~V < 0, 

(3.28) M v+ : = diag (^ +1 , . . . , /O > V > 0, 

(3.29) M p _ := diag ...,//+)< -7/(Re A + |Im A| 2 ) < 0, 

(3.30) M p+ := diag , . . . , //".J > 7?(Re A + |Im A| 2 ) > 0, 
with also 

(3.31) |M p± | = (9(|A|), 



and 

(3.32) M : = 



1 

A a p 



We make a final improvement, following [Go.l-2,MZ.l], by making the rescaling 
lj(x) — > a~ 1 (x)lj(x), rj(x) — > rj(a;)Q;(a:), with a satisfying ODE 

(3.33) a' = —Ijr'jOt. 
This achieves the normalization 

(3.34) Ijr'j = 

while preserving the previously stated properties (3.23)-(3.32). As a consequence, 
we obtain the key property that L R' = 0, i.e., error O e vanishes in the r), z block: 

(3-35) @( v ,z),( v ,z) = 0. 

We could carry out a similar normalization at the Lj, Rj level (with otj now matrix- 
valued) to annihilate all diagonal blocks Qjj; see Lemma 4.9 of [MZ.l]. However, 
we do not require this for our argument. (In fact, we only require vanishing in the 
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lower lefthand corner Z ?7 = of the 77, z block, which is automatic; however, we 
wish to point out early on this more general procedure.) 

First reduction. Observing that coefficients M v _ and M u+ of the "fast trans- 
verse modes" V- and u + are separated by a uniform spectral gap both from each 
other, and from the "slow" modes (p_, 77, z, p+)*, we may apply Proposition 2.6 
twice to reduce (2.13) to decoupled flows on three invariant manifolds, associated 
respectively with z/_, v + , and (p-,rj, z, p+)* , appearing as 

(3.36) v'_ = M v _v_ + 0(<J(e))i/_, 



(3.37) 
and 



v' + = M v+ v + + <D(6(e))v+, 



(3.38) 

where, by (3.25), 
6 e : 



(P~\ 
V 

z 




















1 








A 


a p 





V 








M P+ 



(p-\ 

V 

z 



(p-\ 
v 

z 



is small on the principal block: 



(3.39) 







(77,2), (r),z) 



= 0(S(e)) = 0(e 2 ). 



Observing that the flow is uniformly exponentially decreasing for (3.36) and 
increasing for (3.37), we find that the only possible decaying solution are the trivial 
ones v± = 0, and so we may discard these "fast transverse" equations, leaving us 
with (3.38). Rescaling now by x — > Aex, A — > A/A 2 e 2 , z — > z/Ae, we obtain the 
block-triangular system of equations (1.19)-(1.24) described in the introduction, 
modulo coefficient errors 



(3.40) 
where 
(3.41) 



5 £ £ + e£ e , 



5 £ ( x ) := e- L 5 £ (x/e) = 0(ee~^), 9 > 0, 



with O e uniformly bounded and O(e) on the principal, (r/, z) block), and £ £ = in 
the p± equations and for the (77, z) equations is 



(3.42) 





a p — fj 
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giving 

(3-43) £* ± =(° Q ° {1) y \E*-EU = 0{e- e \*\), 

by Corollary 1.5. Rearranging, we have the error bounds asserted in the introduc- 
tion. Note that the forcing term 

(3.44) ^ A >= • 

„:= e-H'e,,? = 0(e-"W), 

arises in (1.20) as a result of the rescaling z — > z/e. 

Further reductions/Normal forms. Loosely following the introduction, we 
now examine separately each of the three regimes (with respect to the rescaled 
variable A) I. |A| < C, II. C < |A| < Ce' 1 : and III. Ce" 1 < |A| < C^e" 2 , where 
C > is a sufficiently large constant to be determined later. (Recall: In Proposition 
2.9 we have already disposed of the case |A| > C~ 1 e~ 2 .) 

Region I. On region I, the rescaled p±, (r/, z) equations evidently have bounded 
coefficients satisfying condition (ii) of Lemma 2.5, with r/(e) — O(e). Moreover, the 
asymptotic coefficients at plus and minus spatial infinity are block-diagonal, with 
M + and M_ positive and negative definite, respectively, and M = M + O(e), 
where 

M! i) 

at x = ±oo have uniform spectral gap between stable and unstable eigenvalues for 
all Re A > 0; it follows by standard matrix perturbation theory that the stable 
and unstable subspaces for the asymptotic coefficient matrices of the e-approximate 
system are within 0(e) of those of the limiting, e = system, verifying condition 
(ii) of Lemma 2.5, with 77(e) = O(e). Thus, we may conclude from the Lemma 
that the approximate system has a local Evans function D e that is O(e) close to 
the local Evans function (with some chosen basis) D° for the limiting system with 
error terms dropped. As described in the introduction, the limiting p± equations 
decouple, and are readily seen to admit only trivial L 2 solutions, and the limiting 
(r/, z) equation with p± set to zero is just the Burgers eigenvalue equation. Since 
the latter is known to admit no decaying solutions for A = 0, we thus have that 
the limiting Evans function D° has no zeroes on Region I if the Burgers equation is 
strongly stable on Region I, and has a zero of strictly positive real part if the Burgers 
equation is strongly unstable on Region I (possesses a positive real part eigenvalue) . 
By uniform convergence/analy ticity in A, we may conclude the same regarding the 
zeroes of the approximate Evans function D e , and thus also the eigenvalues of the 
original operator about the wave. 
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Region II. On Region II, the eigenvalues 



of the (77, z) ( "Burgers" ) block Mo satisfy 

<-v\M 1/2 <-c 1/2 v<o, 



(3.45) 



/"+>+??|A| 1/2 >C 1/2 ?7>0, 



with also \fip\ 



OdA] 1 / 2 ). Thus, "balancing" the reduced equations by the rescal- 
ing z — > z/IAj 1 / 2 , we obtain well-conditioned basis vectors for this block. At the 
same time, the forcing matrix N of (1.23) becomes now (9(|A| -1 / 2 ) (note: we have 
partially undone the original rescaling z — > z/e), and may be viewed as an error 
term of order C(C -1 / 2 ): uniformly small, but not vanishing as e — > 0. Indeed, 
reviewing the origins of various error terms, we find that this is the dominant er- 
ror term within both of Regions II and III, with error terms in the p± equations 
retaining their originally stated form. 

Diagonalizing M by a well-conditioned change of basis, we thus obtain (in 
rescaled coordinates) an approximately block-diagonalized system 



(3.46) 



\ 



(P 

rj. 
V+ 
P+ 



(P~\ 



= diag(M_,/i-,/i+,M + ) 



V- 



+ 0(|A|- 1 / 2 ) 



(P 

V- 

v+ 

\pj 



where denote the now-separated Burgers modes. Observing that the //^ are 
separated by uniform spectral gap r]\\\ 1 / 2 both from each other and from the C(e|A|) 
blocks M± on all of Region II (indeed, on Region III as well), we may therefore 
apply Lemma 2.6 as second time to further reduce to decoupled equations on three 
invariant manifolds, associated respectively with rj~ , rjp, and (p_,p + ), appearing 
as 



(3.47) 



(3-48) «)'= W + ^C-X- 

and 



(3.49) 




where 5 e is still C(ee" e l x l) as in the original rescaling: note that we are here using 
in a crucial way the fact that the error terms in p± coordinates alone are better than 
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the 0(\X\ x / 2 error terms in the blocks, which were the ones (recorded in the 

uniform description (3.46)) that led to the error/gap ratio lA] -1 / 2 /^] 1 / 2 ~ |A| _1 
determining the tracking angle in Proposition 2.6. 

Noting as usual that 7]p flows are uniformly exponentially growing/decaying, so 
support only trivial decaying solutions, we may discard the first two equations, 
leaving us with the (p_,p+) equation alone. Noting that coefficients M± = O(eX) 
for these "superslow" modes are bounded on Region II, we may then apply Lemma 
2.5 to obtain that suitably chosen local Evans functions for this superslow flow 
converge uniformly (as O(e)) to that for the limiting equations with error terms 
omitted. But, just as in the treatment of Region I, the limiting p± equations decou- 
ple into uniformly exponentially growing/decaying flows, which evidently support 
no nontrivial decaying solutions. Thus, we may conclude, finally, that there exist 
no decaying solutions of the reduced equations, and therefore no eigenvalues for the 
original equations, for A within Region II, and e sufficiently small. 

Region III. The treatment of Region III is similar to but somewhat simpler than 
the treatment of Region II. For, as already pointed out, all steps in the reduction 
to p± equations (3.49) remain valid on this region as well. Noting that, by (1.24), 
M + and M_ have uniform spectral gap of order 6C 2 e on this region, 9 > fixed 
and C > arbitrarily large, while 5 £ remains uniformly O(e), so we may in this 
case apply Lemma 2.6 to obtain the further reduction to decoupled equations on 
two invariant manifolds associated with modes p_ and p + , appearing as 

p'_ = (M_ + 0(C- 2 ^j)p_ 

and 

p' + = [m + + 0(C- 2 5^ p+ )) P+ . 
But, these, having negative and positive definite coefficients 

M± + 0(C- 2 5 e e e p± ) ^ ±9C 2 e + 0{C~ 2 e) 

^ ±6>C 2 e/2, 

are uniformly exponentially decaying and growing, respectively, so support no non- 
trivial decaying solutions. We may thus conclude that there exist no eigenvalues in 
Region III, for e sufficiently small. 

Conclusions. Collecting the results for each regime, we obtain the result of 
Proposition 3.2, and therefore of Theorem 3.1. This completes the analysis of the 
case B = I. □ 

3.2. The general strictly parabolic case. The general case follows quite 
similarly as in the case B = I, in terms of overall argument structure. However, 
there is an important new technical aspect introduced here, in the way that we 
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carry out the linear algebra underlying the initial reduction to "slow" coordinates 
(p±,V,z)- 

For general B, the linearized eigenvalue equations (3.2) become 

(3.50) (B e w')' = (A e w)' + Xw, 
where 

(3.51) B e (x) := B(u e (x)), 
and 

(3.52) A e (x)v := Df(u e (x))v - ^DB(u e (x))v^j (u £ )'(x), 

and the associated integrated equations become 

(3.53) w" = (B e )- 1 A e w' + \(B e )- l w. 
We express these as a family of first-order systems 

(3.54) W = A e (x,X)W, 
where 

(3-55) A ' := (a^V 1 (S £ )-^ £ 

Note that A £ lies within (9(|( , u e ) / |) of dF(u e ), hence is likewise strictly hyperbolic, 
with eigenvalues cij and eigenvectors lj and Tj lying within (9(|( , u <E ) / |) of those of 
dF(u e ). For purposes of our argument, therefore, A e and Df(u e ) are indistinguish- 
able; in particular, cij < —i] < for j < cij > rj > for j > p, and the principle 
eigenvalue a p after rescaling approaches to the Burgers profile f\ with the same rate 
predicted in Proposition 1.5. (This last is the main point, and holds also in the 
more general setting alluded to in Remark 3.4, with (H2) replaced by (H2').) 

For simplicity of exposition, we make the provisional hypothesis: 

(3.56) a(B- l Df(u )) distinct. 



This concerns only the transverse modes j ^ p, since the argument of [MP] shows 
that zero is a simple eigenvalue under assumptions (HI), (H3), and is easily removed 
at the expense of further bookkeeping. The next lemma furnishes the linear algebra 
necessary to carry out the reduction to slow coordinates in the general case; for 
similar arguments, see [ZH,Z.3,MZ.l]. 
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Lemma 3.8. Associated with each of the transverse fields j ^ p of A £ , there exist 
eigenvalues and associated left and right eigenvectors and of A e , analytic 
in A and (A e , B e ) on sufficiently small neighborhoods of A = and (Df(uo), B(uq)), 
with expansions 

(3.57) fi+ = -A/a, + Xtf-ja] + . . . , 



(3.58) L+(x) = -ljB/a 3 ) + . . . , Rf(x) = (j* ) + 
and 

(3.59) nj = 7j - + . . . , 

(3.60) Lj(x) = (0, %/ 7 ,) + . . . , i?-(x) = ) + . 
/or j < p, and expansions 

(3.61) v+ = 7j - + . . . , 



(3.62) L+(x)=(0,V 7j )+... ) i?±(z)= ( 7 ^.)+-- 
and 

(3.63) nj = -A/a, + A 2 /3?/af + . . . 



(3.64) LJ(x)=(l j ,-l j B/a j )+..., Rf(x)=^j+..., 

for j > p. Here, 

7i» • • • » 7 P -i < < < V < 1p+i, ■ ■ ■ , In 

denote eigenvalues of (B e )~ 1 A e , with associated left and right eigenvectors §j, Sj, 
and 13 j := ljB £ rj . 

Dual to the right- and left-invariant subspaces Span j^ p {Rj} and Span j^ p {Lj} 
are two-dimensional left- and right-invariant subspaces spanned by vectors L v , L z 
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and R v , R z that are analytic in X and (A e , B e ) on the same neighborhoods of X = 
and (Df(uo),B(uo)), with expansions 

(3.65) L v = (l p ,l^ + . . ., L z = (0,Sp) + . . . , 
and 

(3.66) «,= (£)+..., *.= (£)+..., 

w/iere 

(3.67) Z,r* = 7j/jS p , Z*s,,- = -IpSj/^fj, for j 
and 

(3.68) Zpr* + Z*s p = 0. 

Proof. By standard matrix perturbation theory, we have analytic dependence 
in (A c ,B e ) and £ of the eigenvalues Xj and eigenvectors lj, fj of the symbol 

(3.69) P e (iO:=-i^ e -£ 2 £ e , 
with expansions 

Aj = -i^oj - £ 2 /3j + . . . , 

(3.70) = /,+..., 

fj = r j + . . . , 

about (D f (uq) , B (uq) , 0, where j3j := ljB e rj. For each j ^ p, aj is bounded from 
zero, and so (3.70) may be inverted, by the Analytic Implicit Function Theorem, 
to yield analytic series for n := i£ and therefore lj, fj in terms of A, (A £ , B e ). This 
verifies expansions (3.57), (3.63) for the "slow" transverse modes; the corresponding 
eigenvector expansions (3.58), (3.64) then follow by the relations 

Lj = (A/V(A + ^.)),^/>7(A + ^), 
= (fjl ^3^3)1 

between lj, fj and Lj, Rj] for similar arguments. 

Expansions (3.59), (3.60) and (3.61), (3.62) for the "fast" transverse modes follow 
in more straightforward fashion by (3.56) and the analytic dependence of isolated 
eigenvalues/eigenvectors of A e at A = 0. Finally, analytic dependence of the comple- 
mentary two-dimensional invariant subspaces follows by duality, whence (see [Kat]) 
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we obtain the existence of analytic bases (L^,L Z ), (R V ,R Z ). By inspection, the 
choices described are valid representatives at A = 0. □ 

Now, similarly as in the case B = I, we make the definitions (3.16)-(3.22), to 
obtain an approximately block-diagonal system satisfying (3.23)-(3.31), but with 
(3.32) now replaced by 

M := L A £ R 

= £o|a=o^ £ -Ro|a=o + (Lo - L |a=o)^[a=0-^0|A=0 
+ L \ X=0 A £ X=0 (R - Ro\\=o) + 0(\X\ 2 ) 

( 3 - 71 ) _ ( XhB- 1 ^ l p s p + XhB^r* + -f p hs p ) 

Xs p B~ 1 r p 7 P + XspB^r* 

0{\X\) 0{\X\)\ 
0(|A|| 7p | + |A| 2 ) 0(|A|)i' 



+ 



where, in the crucial lower lefthand corner of the final (error) matrix, we have used 
the special properties 

^3 72 ^ A |a=o^|a=o = 0, 

£ 2 |A=0 A [a=0 = 7p£z|A=0 

together with the usual 

\Ro - Ro\\=o\, l-^o - ^0|a=oI = 

to obtain the stated bound. The O(-) terms in (3.71) may be expressed more 
precisely as 



(3.73) 



OiA 2 A 
3 A 2 + 4 7 P A O b X 



with coefficients Oj analytic in all parameters. 

Finally, we observe as in [MZ.l] that we may achieve L R' = by a renormal- 
ization L — > a _1 L , -Ro — ¥ Ro&, where a is a 2 x 2 matrix defined by ODE 

(3.74) a' = -L R' a, a| x =o = I- 

At A = 0, LqR'q is initially upper-triangular, hence a is as well. In other words, the 
form of Lq, -Ro are preserved, and (3.74) amounts to the pair of vector normalizations 



together with the condition 



i P K + i*s' p = o 



restricting choices of /*, r*. This does not affect any of the previously stated equa- 
tions or bounds, but only accomplishes again that error © vanish on the rj, z block: 

(3-75) Q (v,z),(v,z) = °- 
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Lemma 3.9. Matrix M defined in (3.71) satisfies 

( o 76) M m _ ( 0{\\\) 0(|A| + e)\ 

(3.76) M -M -^ o{m G(A + e2) J 



and 
(3.77) 

(M - M (±oo)) - (M„ - M o(±0 c)) = ty^JgZ) ) 

for x ^ 0, where 

< 3 ' 78 > *° : =(a//3 a,'//? 

/3 as in the introduction denoting (3 p (uo) = (3 P + 0(e). 
Proof. Taylor expanding about a p = 0, we have 

s P = r p + 1 (a p )a p , 

s p = lpB//3 p + 2 (a p )a p , 
(up to scalar multipliers) and 

1p = a p //3 p + 3 (a p )ap, 

with Oj analytic. From these expansions, together with the fact that A, B, l p , /*, r p , 
r*, s p , s p , and a p , as well as coefficients Oj in (3.73), all approach their asymptotic 
values as x — > ±oo at rate proportional to \u e — u±\ ~ ee _6le ' x ', the result readily 
follows. □ 

From this point on, the argument goes almost exactly as in the case B = I. 
Carrying out the first reduction eliminating fast transverse modes, and rescaling by 
x — > Aex/P, X — > /3A/A 2 e 2 , z — > (3z/Ae, we obtain again the block-triangular system 
(1.19)-(1.24) described in the introduction, with the errors (1.26) there described. 
The latter errors are slightly larger than occurred in case B = I; however, they 
are sufficient that the arguments in each of Regions I III carry through exactly as 
before. This completes the proof of one-dimensional spectral stability in the general 
strictly parabolic case. □ 

Remark 3.10. The crucial estimate of the lower lefthand entry in the final 
(error) matrix of (3.71) is a generalization of the standard observation, in the case 
of a one-parameter bifurcation in A from an isolated semisimple block L(0), R(0) of 
matrix A(A), L(0) and R(0 consisting of genuine eigenvectors with a single common 
eigenvalue /i(0), that 



(3.79) 



LA J R = L(0)A J R(0) + C(A 2 ), 
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which follows from the observation that 

(3.80) (d/d\)LA(0)R(0) + L(0) A(0) (d/d\)R = //(0) (d/dX) (LR) = 0. 

In the case at hand, which is essentially a two-parameter bifurcation in A and a p 
from a Jordan block at A = 0, a p = 0, the computation (3.80) of course breaks 
down in general, but does remain valid for the single, lower lefthand entry of the 
block that corresponds to the pairing of genuine left and right eigenvectors L z and 
R v , respectively, to give the stated result of validity to order (9(|A| 2 + |A||a p |). This 
explains the success of the calculation. Indeed, the entire calculation could have 
been done by expansion about A = 0, a p = 0, giving slightly degraded error terms 
0(\X\ + | a p |) in other entries, since these would in fact have been sufficient for the 
analysis; however, we prefer to carry a p within the computation, since it is no more 
difficult and gives sharper estimates. 

Section 4. One-dimensional stability for relaxation systems. 

Next, we treat small-amplitude stationary profiles of one-dimensional relaxation 
systems 

(4.i) («) =(;, 

\ v Jt \9(^v)J x \q{u,v) 

with q(u, v*(u)) = 0, Re a(q v (u, v*(u))) < as in (1.6)-(1.5), and f(u) = f(u, v*(u)). 
We make the assumptions: 

(HO) /, g, q E C 2 (regularity). 

(HI) a (jyDf, DgY(u, v)*j real, semi-simple with constant multiplicity, and uni- 
formly bounded away from (nonstrict hyperbolicity of relaxation system, plus 
weak subcharacteristic condition). 

(H2) a (Df(u±)) real, simple (strict hyperbolicity of equilibrium system). 

(H3) Re a («±, *>±) + ( °J («±, v±) j < Y+W' 6 > °' f ° r &U 

^eK (linearized stability of constant states). 

(H4) a p = is a simple eigenvalue of Df(uo) with left and right eigenvectors l p 
and r p , and l p D 2 f(r p ,r p ) ^ (genuine nonlinearity for the equilibrium system of 
the principal characteristic field). 

As in the previous section, e > denotes shock strength \u e + — u e _\, and profiles 
(u e ,v £ )(-) are assumed to converge as e — > to (uq,vq) (where necessarily = 
v*(uo)), under which assumptions the center-manifold argument of [MZ.l] verifies 
the assertion of Proposition 1.4, yielding convergence after rescaling of u £ to the 
standard Burgers profile (1.11); see also related analyses of [Yo.Z,FZe]. In the rest 
of this section we establish the following analogs of Theorem 3.1 and Proposition 
3.2. 
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Theorem 4.1. Under assumptions (H0)-(H4), profiles u e are strongly spectrally 
stable (and therefore linearly and nonlinearly orbitally stable [MZ.l]) for e suffi- 
ciently small. 

Proposition 4.2. Under assumptions (H0)-(H4), profiles u e are strongly spec- 
trally stable (and therefore linearly and nonlinearly orbitally stable [MZ.l]) for e 
sufficiently small, if the standard Burgers profile (1.11) is strongly spectrally stable 
for e = 1 and only if the Burgers profile is not strongly spectrally unstable for e = 1 
in the sense that there exists an eigenvalue Re A > 0. 

Theorem 4.1 completes the program of Mascia&Zumbrun [MZ.l], resolving the 
open question of stability of general small-amplitude relaxation profiles. 

Remark 4.3. The "real," or partially parabolic case is closely related to the 
relaxation case, and may be treated by quite similar analysis, following the frame- 
work set up in Appendix A2 of [Z.3]; see [Z.6]. In the real viscous case, satisfactory 
small-amplitude stability results have already been obtained by energy methods 
[HuZ]. 

4.1. The 2x2 case. For clarity, we first carry out the 2x2 case /, g, q, f e M 1 , 
which was treated by energy methods in [L.3]. The linearized eigenvalue equations 
about profile (u e ,v e ) are 

(4.2) H:))'^-^:). 

where 

(4.3) A £ (x):=(( u lA(u',v e )(x), Q"(x):=(° 

\9u 9v J \Qu 

Following [Z.3], we make the (nonsingular, by (HI)) change of variables 

(O-'C) 

to obtain a first-order system 

= (Q e -A/)(^)- 1 (Q 

in nondivergence form. 

So far, this is analogous to the "flux transform" of [Go. 2]. If we were now to 
convert to a single, second-order equation for /, then the analogy would be exact, 
since the equations would then clearly not support any L 2 (decaying) solutions at 



J (u^v^)(x). 



(4.4) 
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A = 0. Instead, we take a different approach, keeping the convenient variables /, g, 
and rescaling as 

(4-6) /-/A, ~g-h 

to obtain 

(4.7) (J^j = diag(A- 1 ,l)(Q e -A/)(^)- 1 diag(A,l) 

in place of (4.5). Since /' = Xu in the original equations (4.2), we see that the 
rescaled / variable is, in the case of decaying solutions, given by f(x) = u{z)dz, 
i.e., is exactly the integrated variable of standard energy analyses [L.3,Liu,HuZ]; in 
particular, the rescaled equations (4.7) have the favorable property that they do 
not support L 2 (decaying) solutions at A = 0. 

Remark 4.4. Equations (4.7) are neither pure flux nor pure integrated form, but 
a new form intermediate to the two, which we might be called "balanced flux form." 
This form will prove useful again in the multidimensional case treated in [PZ]. (All 
three forms (flux, integrated, and balanced flux) coincide in the one-dimensional 
viscous case, so this distinction is not apparent, or necessary, here.) 

By direct calculation, we have 



A' 1 = (det A)' 1 
QA' 1 = (det A)' 1 ^ 



9v fv 
9u fu 





Qu(jv Qv(ju Qufv ~\~ Qvfit 



q v ( 
det A V ~9u «i 
where g(u) := g(u, v*(u)). Combining, we find that 

^ \ ( 9v -f v 

~9u fu 

\f v /q v \ + fO(\X\) 



(Q-X)A-i = ( d etAr( qv (_° gu a ° i )-A( 



det^V"^ oi J 0(\X\) 

and therefore 

diag (A" 1 , 1) (Q - X)A~' diag (A, 1) = 

Qy ( f v /q v \, f 0(\X\) 
det A\-g u X ai J + \0(\X\ 2 ) 0{\. 

in (4.7). 
Let 

(4-8) /?i := -f v q~ 1 (g u - q^Qufu) = -fvQ^teu - 9« 1 9uOi) 

denote the effective, "Chapman-Enskog viscosity" defined in (1.15). A consequence 
of (H3) is that (3 > 0; see, e.g. [Ze], or Appendix A.2 of [MZ.l]. 
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Lemma 4.5. There hold det A = —f v gu + g v a± = —f v g u + 0{t) and 

Pi = -qv 1 fvgu + fvqv 2 Quai = -g^det A + {q~ X g v + A,<?~VW 

= -g-Met A + 0(e). 

In particular, det A < and sgn f v = sgn g u ^ 0, for e sufficiently small. 

( f f \ 

Proof. A column operation reduces det A to det ~ J > from which the 



g u g v _ 

first assertion follows. Combining with (4.8), we obtain the second assertion. □ 

Remark 4.6. Condition det A < for a\ near zero may be recognized as the 
subcharacteristic condition, which in the 2x2 case is equivalent to (H3). Condition 
f v 7^ is the "genuine coupling" assumption of Liu [L.3]; as we have just shown, it 
is in fact a consequence of (H3), and not a separate assumption. 

Now, setting 

:= P 1 A/q v (0), 

rescaling x — > Aex/(3, A — > /3A/A 2 e 2 , and balancing / — > —g u (3f/Ae, we obtain the 
normal form 



g J \X ai J \g J \g 

where, by Lemma 4.5, error term $ satisfies bounds (1.22)-(1.26) of the introduc- 
tion: that is, just an approximate Burgers eigenvalue equation. From this form we 
may immediately conclude the stability results, by a simplified (because transverse, 
p terms are absent) version of the arguments of the previous section. This completes 
the analysis of the 2x2 case. □ 

4.2. The general case. Now, consider the case of general n, r > 1, starting 
as in the previous subsection with the linearized eigenvalue equation written in the 
"balanced flux variables" introduced in the previous subsection: 

(4.9) W = A e W, 



A e := diag(A- 1 / n ,/ r )(Q e -A/)(^ e )- 1 diag(A/ n ,/ r ) 
( 4 - 10 ) _(ln OA i/AJn 

q u q v ) V o i r 

A e , Q e as defined in (4.3). 
Define 



(4.11) (H,H):=(q u ,q v )A- L , 
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H e R rXn , H e M rxr . As discussed in the treatment of existence given in Appendix 
A.l of [MZ.l], the r x r matrix H has a special role as the coefficient matrix of 
the linearized traveling wave ODE expressed in flux coordinates. For simplicity of 
exposition, similarly as in the viscous case, we make the provisional hypothesis: 

(Al) aH(uo) distinct; 

again, this is easily removed at the expense of further bookkeeping. A consequence 
of dissipativity, (H3), is that on a sufficiently small neighborhood of the base point 
wo, H(u) has a single eigenvalue 

(4.12) lk = 0(e), 

with the (r — 1) remaining eigenvalues uniformly bounded from zero as e — > 0: 

(4.13) 7i,...,7 fc _i < -r) < < T) < 7fc+i,...,7 r ; 

we denote the associated left and right eigenvectors by §j and Sj, respectively. 
At A = 0, A e reduces to 

(4.H) (» * 

where 

(4.15) (E,E):=(^-I n ,0yAT 1 - 
More generally, we have 

(4.16) A«(x,A)=( A // A ^ H * xp ), 

where E, E, H, H are as defined in (4.15) and (4.11), and 

(4.17) (F,F):= (O,-/,)^)- 1 . 

Lemma 4.7. Associated with each eigenvalue jj, j ^ k of H , there exist "fast" 
eigenvalues n - and associated left and right eigenvectors L J - and R - of A e , analytic 
in X and (A e , Q e ) on sufficiently small neighborhoods of A = and (A e , Q c )(uq, vq), 
with expansions 

(4.18) ^ = 7j + ..., 
and 



(4.19) if =(0, %)+..., If {■■'f' / ) 



+ ..., 
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where eigenvalues and associated left and right eigenvectors of H, 

and E is as defined in (4-15). 

Likewise, associated with each of the transverse fields j 7^ p of A £ , there ex- 
ist "slow" eigenvalues (ij and associated left and right eigenvectors L s - and R s , of 
A e , analytic in A and (A e , Q e ) on sufficiently small neighborhoods of A = and 
(A e , Q e )(uo,vo), of A = and (A, Q)(uq,vq)), with expansions 



(4.20) ^ = -A/a, + X 2 pj/a 3 j + . . . , 

(4.21) /-•(•'•) (/,./,) • .... RZ(x)=( r A+..., 

where (3j := ljB*rj, with £?* the effective, " Chapman- Enskog" viscosity defined in 
(1.15) of the introduction, and tj satisfies 

(4.22) ljE + tjH = 0. 



Dual to the right- and left-invariant subspaces Span j^k{Rj} © Span j^ p {Rj} 
and Span j^k{Lj} © Span j^ p {Lj} are two-dimensional left- and right-invariant 
subspaces spanned by vectors L v , L z and R v , R z analytic in A and (A e , Q e ) on the 
same neighborhoods of A = and (A e , Q e )(uo, vq), with expansions 

(4.23) L v = (l p , /*)+..., L z = (o, s p ) + . . . , 
and 

(4.24) B „=( r *)+..., fl,= (^)+..., 
where 

(4.25) Ijr* = —tjSp, l*Sj = —i^lpEsj, forj^p 
and 

(4.26) Z p r* + hs p = 0. 

Proof. This follows similarly as in the proof of Lemma 3.8, namely, analytic 
variation of the fast transverse modes follows by their spectral separation, while 
analytic variation of slow transverse modes follows by inversion (and appropriate 
conjugation) of the associated dispersion relations 

Xj = -iiaj - i 2 l3j + . . . , 

(4.27) C, 

n j =n j + ..., 
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of the symbol P e {i^) := Q £ — i^A £ . The calculation of expansions (4.27) may be 
found, e.g., in Appendix A. 2 of [Z.3] or Appendix A. 2 of [MZ.l]. Analytic variation 
of the two-dimensional complementary subspaces again follows by duality. □ 

As in the treatment of the viscous case, we make the definitions (3.16)-(3.22), to 
obtain an approximately block-diagonal system satisfying (3.23)-(3.31), with (3.32) 
replaced by the following analog of (3.71): 

M := L A £ R 

= L |A=0-^ £ -Ro|A=0 + (Lo - £o|A=o)^fA=0-Ro|A=0 
+ L \x=oM x=0 (R - -RoiA=o) + C(l^| 2 ) 



(4.28) 



+ 



I A=0 jf ^[A=0^ xx - -rt0|A=0j 

0(A) l p Es p + -f p hs p + 0(X)\ 
Xs p Hr p 7p + 0(A) J 

( 0(\\\) 0{\\\)\ 
x O(|A|| 7p | + |A| 2 ) 0(\X\) J ' 



where all O(-) terms are analytic in A, A, Q. Here, as in the viscous case, we 
have used (3.72) to improve the crucial bound in the lower lefthand corner of the 
final, error matrix; more precisely, we have description (3.73) for this error bound, 
with coefficients Oj analytic in all parameters. As in the viscous case, we may 
achieve LoR' = by a renormalization leaving the form of the equations otherwise 
unchanged; thus, error O vanishes on the ?y, z block: 

(4-29) ©<.m),Om) = 0- 

Lemma 4.8. Matrix M defined in (4-28) satisfies 



(4.30) Mo-Mo=(^' c ) 



0{\\\) 0(|A| + e) 
0(X + e 2 ) 



and 
(4-31) 

(Mo - M (±oo)) - (Mo - M (±oo)) = ( ^[^I^l j ^ + 



for x ^ 0, where 

(4.32) M 



1 

A//3 a p /(3 



(3 := (3 P (0), and as usual (3j := ljB*rj with £?* as defined in (1.15). 

Proof. Taylor expanding about a p = 0, we have up to scalar multipliers (see 
Appendix A.l, [MZ.l], Claim just below equation (8.7)): 

(4.33) s p = g u r p + Oi(\a p \)a p , 
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(4.34) s p = -l p f v q v 1 //3 p + 2 (\a p \)a p , 

and (by a standard matrix perturbation calculation for a simple eigenvector: see 
(8.14)-(8.15) of Appendix A.l, [MZ.l]). 

(4.35) 7 P = a p /(3 p + 3 (\a p \)al, 

where Oj are analytic. 

Furthermore, at a p = 0, there hold the relations: 

(4.36) (e,e)^=t p , 



(4-37) (P,f) ^=- q -\ u r p 

(see (8.8), Appendix A.l, [MZ.l]), and 

(4.38) s p (h,H^ = (-J p /A0) 

(see (8.10), Appendix A.l, [MZ.l]), from which we obtain the key facts: 

(4.39) Esp = r p , Fs p = -q~ 1 q u r p , s p H = -l p //3. 

From expansions (4.33)-(4.35), the relations (4.39), and the fact that A, /3 P , Z p , 
/*, r p , r*, Sp, Sp, and a p , as well as coefficients Oj in (3.73), all approach their 
asymptotic values as x — > ±oo at rate proportional to \u e — u±\ ~ ee _6,e l x l), the 
result readily follows. □ 

Remark. Note that the natural choice of basis vectors L z and R z given in 
Lemma 4.7 automatically accomplishes the desired rescaling to the viscous case, so 
that we need not do this "by hand" as in the treatment of the 2x2 case above. 

From this point, the argument goes through exactly as in the general viscous 
case. Carrying out the first reduction eliminating the r — 1 fast transverse modes, 
and rescaling by x — > Aex/fi, A — * /3A/A 2 e 2 , z — > /3z/Ae, we obtain as in the general 
viscous case the block-triangular system (1.19)-(1.24) described in the introduction, 
with errors (1.26), and so the remaining arguments go through as before. This 
completes the proof of one-dimensional spectral stability in the general relaxation 
case. □ 



Section 5. Multidimensional stability for parabolic systems. 
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Finally, we briefly discuss the extension of the above analysis to the multidi- 
mensional viscous case; details will be given in [PZ]. Consider a sequence of planar 
stationary profiles u e (xi) of a strictly parabolic viscous conservation law of form 

(5-1) u t + J2f j (u) Xj = Y,(Bi k (u)u Xk ) Xj , 

j jk 

lying in a neighborhood U of a particular state u , where x G R d , u, / J G M n , and 
B jh G R nXn . We make the assumptions: 

(HO) p,B^ k G C 2 (regularity). 

(HI) Re a(J2B jk (u)^ k ) > for f G R d (strict parabolicity) . 

(H2) There exists A (•), symmetric positive definite and smoothly depending on 
u, such that A G Dfi(u) is symmetric for all 1 < j • < d. (simultaneous symmetriz- 
ability, =>- nonstrict hyperbolicity) . 

(H3) Re a(-iZDP(u ± )Cj -E^(«±)^) < -#l£| 2 , £ e R d , for some 9 > 
(linearized stability of constant states). 

(H4) (i) a p = is a simple eigenvalue of L>/ 1 (w ) with left and right eigenvectors 
l p and r p , and l p D 2 f 1 (r p ,r p ) ^ (genuine nonlinearity of the principal character- 
istic field in the normal spatial direction xi). 

(ii) r p is not an eigenvector of any Df^ := X^yi ^j-^/ J f° r £ e R d_1 ^ 0, and, 
in the intermediate case 1 < p < n: 

(5.2) (r p , A°A^fl{A 1 )- 1 ILA^r p ) ^ 

for all | G ^ 0, where A^ := Y.j^ P D f j iji Al '■= Df 1 , and n denotes the 

projection complementary to the eigenprojection of A 1 onto r p ( "nonresonance," or 
genuine coupling in the sense of [Me.l-4,FZ,Z.5-6]). 

Here, e > as usual denotes shock strength \u e + — itf_|, and profiles u e (-) are 
assumed to converge as e — > to -u - Since the profile problem is restricted to the 
normal direction, the one-dimensional analysis of Majda&Pego [MP] again yields 
Proposition 1.4, with convergence after rescaling of u e to the standard Burgers 
profile (1.11). Linearizing about u e , we obtain a linearized evolution system (1.8) 
with coefficients depending only on xi; taking the Fourier transform in directions 
x := (2:2, • • • , Xd), we then obtain for each fixed e a family of one-dimensional equa- 
tions 

(5.3) v t = L^v, 

indexed by f G lR d_1 . 

Definition 5.1. Following [Z.3], we define strong spectral stability (condition 
(DP) of the reference) in the multidimensional case as 



(5.4) 



Re a p (L c ) < Ofor (f, A) ^ (0,0) 
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equivalently, the generalized eigenvalue equation (L| — X)w = admits no L 2 solu- 
tions for Re A > and (£, A) ^ (0, 0). 

Under hypotheses (H0)-(H4), the results of [Z.3] show that strong spectral spec- 
tral stability implies linearized and nonlinear stability with sharp rates of decay (as 
discussed in [PZ], our conditions (H0)-(H4) imply conditions (H0)-(H7) assumed 
in the reference). In [PZ] we establish: 

Theorem 5.2. Under assumptions (H0)-(H4), profiles u e are strongly spectrally 
stable (and therefore linearly and nonlinearly orbitally stable [ZH,Z.2]) for e suffi- 
ciently small. 

Remark 5.3. As in the inviscid small amplitude case [Me.l-5,FM,Z.5-6,FZ], the 
condition of simultaneous symmetrizability (H2) may be relaxed under appropriate 
alternative structural conditions; however, so far there has been identified no such 
conditions of simple and general application. Likewise, the nonresonance condition 
(H4)(ii) might possibly be relaxed; however, the degenerate case that this condition 
fails would appear to require a considerably more detailed analysis than the one 
carried out here. In any case, these conditions are satisfied for most systems arising 
in physical applications; see [Me. 5, Z. 3, Z. 5-6]. In the extreme shock case, (H4)(ii) is 
equivalent to a simpler and easily verified condition identified by Metivier [Me. 1-5]: 
that the principal eigenvalue a p of the full symbol Y^j=i ^j^f^ De a strictly convex, 
resp. concave, function of £ = (£i, • • • , £d), according as p = 1 or n; see [FZ] or 
(implicitly) [Me.l]. Note that (5.2) is automatically satisfied in the extreme shock 
case. The case of "real," or partially parabolic, viscosity, or relaxation may be 
treated similarly, following the model of the present, one-dimensional analysis. 

Similarly as in the one-dimensional case, Theorem 5.2 is established by showing 
that stability is equivalent, modulo pure imaginary eigenvalues as in Theorems 3.2 
and 4.2, to stability of an appropriate canonical system. However, this system is 
not a 2 x 2 system as in the scalar multidimensional case, but a 3 x 3 system 
corresponding to the full slow flow of an appropriate 2x2 viscous conservation law, 
namely, the coupled Burgers- linear degenerate equations, 

^ (:) t + ( !, 2 2 ) + (:) =i>-(:) - 

considered with the family of profiles 



where a > and £? jfc are constant, diagonal, and satisfy the uniform parabolicity 
condition 



(5.6) 




(5.7) 



jk 
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for all f := 6) G ^ 2 , with B\\ = 1. 

That is, m multiple dimensions, small amplitude behavior is not captured by any 
scalar model, but generically requires a 2 x 2 model for its expression. This fun- 
damental observation was made in the inviscid setting by Metivier [Me.l-5,FM] in 
his study of the weak shock limit, and is due to the new phenomenon of "glancing 
modes" arising in the multi-dimensional case; see also [Z.5-6,FZ] for related dis- 
cussions in the inviscid case. Accordingly, the analysis of the reduced equations 
on the "slow-superslow" manifold becomes considerably more complicated in mul- 
tiple dimensions; in particular, we can no longer rely on simple Sturm-Liouville, 
or maximum principles to conclude stability of the reduced system. The following 
proposition, though it does not (because the equations (5.5) are not scale-invariant) 
directly imply stability of the reduced system, nonetheless indicates at an intuitive 
level why the analysis can still be carried out, to yield a full stability result and 
not only a reduction. Note that the reduction to constant coefficients a, A 2 , B^ k 
(achieved as in the one-dimensional analysis by an application of Lemma 2.5) is 
crucial in the proof. 

Proposition 5.4. Shock profiles (5.6) of any amplitude e are strongly spectrally 
stable as solutions of (5.5), assuming only the partial parabolicity condition 

(5-8) Y,^ Bik ^ e \tf(l o)' 9>0 - 

jk v J 



Proof. We use a straightforward energy estimate modifying that used in [MN,Go.l] 
to treat the one-dimensional genuinely nonlinear scalar case. Linearizing (5.5) about 
a stationary profile (5.6), we obtain the generalized eigenvalue equation 

(5.9) (A + i&A 2 + ilB 22 )w + (Aw)' + i&{B 12 + B 21 )w' = B 1 V ', 
where 

(5.10) *:=(V°J 

here, as usual, "' denotes d/dxi. 
Examining the symbols 



of the constant coefficient limiting equations of (5.9) as x\ — > ±oo, we find from 
(5.8) that 



(5.11) Re P± > 9\Z\ : 



1 
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whence, for £ 7^ 0, P± has no pure imaginary eigenvalues A unless associated with an 
eigenvector of form (0, 1)*, in which case we find by direct calculation that £2 = and 
A = — a£i and (5.9) decouples into the one-dimensional Burgers eigenvalue equation 
and a trivial, constant-coefficient equation. In either case, we find from standard 
considerations (see, e.g., [He]) that w and all derivatives decay exponentially at 
spatial infinity for w G L 2 and (£, A) 7^ (0, 0). 
Similarly, from 

Re (X + i^A 2 +e 2 B 22 ) = Re (A)/ + £ 2 2 Re (B 22 ) > 2 £ 2 2 (J °) 

for B 2 > 0, by (5.7), we easily find that (A + i&A 2 + i\B 22 ) is invertible for all 
real £ 2 and Re A > such that (£ 2 , A) 7^ (0, 0), whence w(z)dz = 0; for, this 
implies that the kernel, if it exists, is spanned by (0, 1)*), yielding £2 = A = by 
direct calculation. Combining these observations, and integrating (5.9) from —00 
to +00, we obtain 

/+00 
w(z)dz = 
-00 

and therefore w(z)dz = 0. 

Thus, we may introduce the integrated variable W(x\) := J^ 1 oo w(z)dz, with 
W E L 2 satisfying the integrated eigenvalue equation 

(5.13) (A + ii 2 A 2 + i 2 2 B 22 )W + AW + ii 2 {B 12 + B 21 )W = B X1 W" ', 

with W and derivatives again decaying exponentially at spatial infinity. Taking the 
real part of the complex inner product of W against (5.13), integrating by parts, 
and using (5.7), we obtain the estimate 

Re X\W\ 2 - J u'\Wl\ 2 dx = -{W, B X1 W) + (i&W, B 12 W) 

+ (W, B 12 i&W) + {i&W, B 22 it 2 W) 
<-0{\W[\ 2 + ii\Wt\ 2 ), 

6 > 0, a contradiction of the fact that u' < 0, unless \W[\ = ^l^il = 0, hence 
\W\\ = 0, and Re A = 0. But, this case may be disposed of by observing that (5.9) 
then reduces to a single equation in W 2 , which is constant-coefficient and therefore 
supports no L 2 solutions. □ 
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